Spatial and Spatio-temporal Epidemiology

Introduction Understanding provider preferences for the design of HIV treatment packages could enhance the implementation of programs to support the adoption of long-acting injectable antiretroviral therapy (LAI ART) by people living with HIV who are interested in initiating this treatment modality. Methods We recruited providers from New York City (NYC), Rockland, Putman, and Westchester County Ryan White Part A Medical Case Management (MCM) programs to complete a discrete choice experiment (DCE) containing twelve tasks with two alternatives and an opt-out option, with additional survey questions about implementation readiness and choice motivations. The alternatives included four attributes--Type of ART Medication (monthly or bimonthly LAI ART), Service Location and Mode, Support for Clients, and Rewards for Clients--with 2-4 levels each. We ran latent class multinomial logit analyses (LCA) with 1-5 classes to estimate preferences and explore hypothesis-free preference heterogeneity. We estimated attribute influence using relative importances and preferences using zero-centered part-worth utilities for each level. Results One hundred seventy-seven providers completed the survey (July 2022-January 2023). About half (52%) were 40-59 years old, 72% identified as women, and the plurality (41%) identified as Latino/a. We chose the two-group LCA solution. Bimonthly LAI ART was preferred over monthly LAI ART overall and in both groups. Group 1 (n=45) preferred more traditional adherence supports (e.g., injections at the clinic by appointment, injection appointment reminders) whereas Group 2 (n=132) preferred more client-centered supports (e.g., injections at home by appointment, free transportation to injection appointments if at a clinic). Both groups preferred higher monetary value gift cards for clients for every on-time injection. The top-ranking motivations indicated that participants prioritized patient convenience over job satisfaction and administrative or financial feasibility for the agency. The scores for all implementation measures indicate readiness to implement LAI ART in both groups. Conclusions Our implementation science-focused study suggests that providers of MCM services in NYC and surrounding counties are motivated to offer services to support clients' access and adherence to LAI ART. More work is needed to understand how programs have, in fact, integrated supports for LAI ART into their services.

BackgroundDrowning remains a major global public health challenge. This study examined whether the timing and trajectories of urbanisation--beyond the current built environment--are associated with subnational drowning mortality. MethodsWe linked satellite-derived measures of built-environment change (GHSL), population crowding (WorldPop), surface water exposure (JRC Global Surface Water), and infrastructure proxies (VIIRS/DMSP nighttime lights) to GBD 2021 drowning mortality estimates across 203 ADM1 regions in 12 countries (2006-2021; 3,248 region-year observations). Temporal predictors captured recent expansion, development "newness" ([≤]10-year built share), acceleration/volatility, and a crowdingxgrowth interaction. We screened predictors using LASSO (10-fold cross-validation) and fitted mixed-effects models with region random intercepts. Distributed-lag models tested temporal precedence and development age, and income-stratified models assessed heterogeneity. ResultsAdding temporal predictors improved fit beyond contemporaneous built-environment measures ({Delta}AIC=177; {Delta}BIC=147). In adjusted models, crowdingxgrowth was strongly positively associated with drowning mortality, and a higher share of recent development was associated with higher mortality. Lag models showed a development age gradient: older built environment was most protective. Associations differed by income group, with several key coefficients reversing sign across strata. DiscussionDrowning mortality appears shaped by development histories as well as present-day conditions, with risk concentrated in rapidly changing, dense settings and the newest built environments. Cross-context heterogeneity suggests mechanisms and prevention priorities are unlikely to be uniform. ConclusionsDevelopment timing and trajectories help explain subnational drowning mortality beyond current built form alone. Prevention and planning should prioritise transition-period safety strategies in newly developing and rapidly densifying areas.

Drowning remains a major global public health challenge, yet how built environment characteristics shape population-level drowning risk remains poorly understood. This study linked satellite-derived built environment data to subnational drowning mortality estimates across 203 regions in 12 countries from 2006-2021. It found that built environment associations with drowning mortality are complex, non-linear, and shaped by development context. Urban extent was strongly protective, while built area near water showed protection overall but increased risk when combined with high population crowding. Almost all drowning mortality variance occurred between regions rather than within regions over time, indicating risk is predominantly determined by place-based characteristics. Income-stratified analyses revealed profound heterogeneity: crowding was protective in low-to middle-income settings but near-null in high-income regions, while waterfront development captured very different realities across contexts. These findings highlight the importance of tailoring drowning prevention strategies to local built environment configurations and development contexts.

ObjectivesRising temperatures are a major climate-related hazard for U.S. workers, increasing heat-related illness and a broad range of occupational injuries through indirect pathways often overlooked in economic evaluations. We examined the association between temperature and occupational injury and illness and quantified heat-attributable injuries (including illnesses) and costs in New York State. MethodsWe conducted a time-stratified case-crossover study of 591,257 workers compensation (WC) claims during the warm season (2016-2024). Daily maximum temperature was linked to injury date and county and modeled using natural cubic splines, with effect modification by industry and worker characteristics. ResultsInjury risk increased with temperature, becoming statistically significant at approximately 78{degrees}F. Relative to 65{degrees}F, injury odds increased to 1.06 (95% CI: 1.01-1.10) at 80{degrees}F, 1.12 (1.07-1.18) at 90{degrees}F, and 1.17 (1.11-1.23) at 95{degrees}F. Overall, 5.0% of claims (2,322 annually) were attributable to heat. At temperatures [≥]80{degrees}F, an estimated 1,729 excess injuries occurred annually, generating approximately $46 million in WC costs. An estimated $3.2 million to $36.1 million in medical expenditures were associated with incomplete claims, likely borne outside the WC system. ConclusionsThese findings demonstrate substantial economic costs not fully captured within WC and support workplace heat protections as a cost-containment strategy that can reduce health care spending and strengthen workforce resilience.

IntroductionDrowning risk begins with water exposure, yet population-water relationships have rarely been quantified at scale using environmental measures. This study explored whether satellite-derived data was associated with subnational drowning mortality and whether associations differed by country income level. MethodsWe linked Global Burden of Disease (GBD 2021) age-standardised drowning mortality rates to satellite-derived exposures for 212 subnational regions across 12 countries (2006-2021; 3,392 region-years). Exposures were extracted via Google Earth Engine and standardised. Gamma-log generalised linear mixed models included region random intercepts and year fixed effects. Income-stratified models were estimated separately. Supplementary models assessed maritime vessel activity. ResultsNear-water population percentage was the strongest correlate of drowning (IRR 1.40; 95% CI 1.33-1.47). Permanent water coverage was protective (IRR 0.80; 0.73-0.88), as were nighttime lights (IRR 0.96; 0.95-0.97) and hot days [≥]30{degrees}C (IRR 0.95; 0.92-0.99). Mean temperature (IRR 1.17; 1.11-1.23) and precipitation (IRR 1.03; 1.01-1.04) were positively associated. Near-water effects were consistent across income strata (LIC 1.25; MIC 1.31; HIC 1.24), while other predictors showed weak or inconsistent within-strata associations. Vessel activity was modestly associated with drowning in Global Fishing Watch models (IRR 1.05; 1.01-1.09) but not in Synthetic Aperture Radar models. DiscussionSatellite-derived indicators can characterise drowning risk at scale, with population proximity to water emerging as a robust cross-context correlate. Protective associations for permanent water suggest landscape configuration may shape risk beyond proximity alone, highlighting geospatial datas value for targeting prevention where surveillance is limited.

Socioeconomic inequalities in health and well-being are a major public health concern, particularly in ageing populations. Education is a key determinant shaping multiple aspects of health outcomes. We used cross-sectional data from wave 9 of the German sample (n=4,148) of the Survey of Health, Ageing and Retirement in Europe (SHARE) to test whether formal education is associated with well-being in later adulthood, with health literacy, self-rated health, and preventive health behaviours as possible mediators. Our results showed that education was positively associated with greater well-being, but only via indirect pathways. Specifically, self-rated health, health literacy, and fruit and vegetable consumption mediated the relationship between education and well-being accounting for 54.7, 24.7, and 12.6 percent of the total effect, respectively. In addition, there were significant positive correlations between education and health literacy, as well as high-intensity physical activity, daily fruit and vegetable consumption, more preventive health check-ups, and less smoking. In contrast, alcohol consumption was more common among those with higher levels of education. All health behaviours and health literacy were correlated directly or indirectly (i.e., mediated by health) with well-being. These findings highlight the importance of examining indirect pathways linking education to well-being in later life. Interventions aimed at improving health literacy and promoting healthy behaviours may help reduce educational inequalities in quality of life among older adults.

Drug overdose deaths in the United States reached record levels during the fentanyl era before recently declining. A plausible hypothesis is that a sudden drop in fentanyl purity beginning in 2023 caused the downturn in overdose mortality. We evaluated this hypothesis by replicating a published analysis with regional overdose data, using models that account for time trends and autocorrelation, and negative control indicators to test for spurious correlation. When fentanyl purity was rising, the national purity series did not track overdose increases in most regions and showed only a modest association in the West. When both purity and mortality later declined, the observed associations were also seen with unrelated macroeconomic indicators that shared the same time pattern. National fentanyl purity alone does not provide a sufficient explanation for recent overdose declines.

Understanding the dynamics of HIV epidemics is important to control them effectively. Classical methods that mainly rely on occurrence data are limited by the fact that an unknown part of the epidemic eludes sampling. Since the early 2000s, phylodynamic methods have enabled the estimation of key epidemiological parameters from virus genetic sequence data. These methods have the advantage of being less sensitive to partial sampling and to provide insights about epidemic history that even predates the first samples. In this study, we analysed 2,205 HIV sequences from the French ANRS PRIMO C06 cohort. We identified and were able to reconstruct the temporal dynamics of two large clades that represent the HIV-1 epidemics in the country. Using Bayesian phylodynamic inference models, we found that the first clade, from subtype B, originated in the end of 1970s, grew rapidly during the 80s before decreasing from 2000 to 2015 and stagnating since then. The second clade, from circulating recombinant form CRF02_AG, emerged and spread in the 80s, grew again in the early 2000s, before declining slightly. We also estimated key epidemiological parameters associated with each clade. Finally, using numerical simulations, we investigated prospective scenarios and assessed the possibility to meet the 2030 UNAIDS targets. This is one of the rare studies to analyse the HIV epidemic in France using molecular epidemiology methods. It highlights the value of routine HIV sequence data for studying past epidemic trends or designing public health policies.

Tick-borne diseases are now reported from nearly every county in Illinois, and three vector tick species (Amblyomma americanum, Dermacentor variabilis, and Ixodes scapularis) are of particular concern because these are responsible for most of the tick-borne disease transmission in the state. However, active surveillance is patchy, many counties have little or no sampling, and there is no statewide, quantitative map of relative abundance that can be used to anticipate risk in unsampled areas. To address these gaps, we developed Bayesian hierarchical spatial models to estimate the county-level abundance of these three vector tick species in Illinois. Using active surveillance data from 2019-2022, we modeled county-level abundance as a function of climate, land cover, forest fragmentation, and deer habitat suitability. Spatial dependence was captured using a Besag-York-Mollie 2 (BYM2) prior implemented in INLA, along with spatial 5-fold cross-validation to assess predictive performance. A. americanum showed the highest predicted abundance in southern and central Illinois, D. variabilis was widespread but diffuse, and I. scapularis was concentrated in northern and selected central counties. Together, these models provide the first spatial, statewide, uncertainty-aware assessment of tick abundance in Illinois, highlighting priority counties where surveillance lags disease risk.

Objective: To describe the clinical characteristics of term neonates with neonatal bacterial meningitis (NBM) and explore the association between different pathogens and imaging complications, providing clinical evidence for early identification and individualized management. Methods: A retrospective study was conducted on 531 term neonates diagnosed with NBM admitted to the Capital Institute of Pediatrics from 2013 to 2025. Demographics, clinical manifestations, laboratory parameters, etiological results, imaging complications and treatment measures were collected. Patients were divided into favorable/adverse discharge outcome groups and pathogen-positive/negative groups. Statistical analyses were performed using appropriate tests, and Cramers V coefficient was used to analyze the association between pathogens and imaging complications. Results: (1) The most common clinical manifestations were abnormal body temperature (79.85%), altered consciousness (55.18%) and jaundice (46.52%). CSF/blood culture was positive in 133 cases (25.05%), with Escherichia coli (27.07%), group B streptococcus (17.29%) and Staphylococcus species (16.54%) as predominant pathogens. The overall incidence of imaging complications was 22.22%, mainly hydrocephalus (5.84%), subdural effusion (4.90%) and encephalomalacia (2.64%). (2) Adverse discharge outcomes occurred in 107 cases (20.15%). Compared with the favorable group, the adverse group had higher incidences of convulsions, altered consciousness, anterior fontanelle bulging, abnormal muscle tone and primitive reflexes (all P<0.001), more obvious laboratory abnormalities (higher CRP, CSF leukocytes and protein, lower CSF glucose, all P<0.05), higher culture positive rates and greater need for adjuvant therapy (all P<0.001). (3) Pathogen-positive patients had higher imaging complication rates. Gram-negative infections were associated with higher hydrocephalus and subdural effusion rates, while Gram-positive infections had higher brain abscess risk. Specifically, Escherichia coli correlated with hydrocephalus and subdural effusion; group B streptococcus with cerebral infarction and encephalomalacia; LM with intracranial hemorrhage and brain abscess; negative cultures correlated with no imaging complications (all P<0.05). Conclusion: Term NBM neonates have non-specific manifestations, mainly abnormal body temperature and altered consciousness. Predominant pathogens are Escherichia coli, group B streptococcus and Staphylococcus species, with hydrocephalus and subdural effusion as common imaging complications. Adverse outcomes are associated with severe symptoms, obvious laboratory abnormalities and higher pathogen positivity. Specific pathogens correlate with distinct imaging complications.

PurposeTo characterize maternal transport patterns in Iowa, a state with levels of maternal care and without formal perinatal regions, and assess whether transport decisions reflect efficient, risk-appropriate coordination. MethodsWe analyzed 2010-2023 Iowa birth records, which included 2,251 maternal transports between obstetric facilities across 106 unique routes. We characterized transport patterns and applied a community detection algorithm to identify "communities" of obstetric facilities that disproportionately transport among themselves. FindingsSuburban and rural counties have elevated transport rates compared to urban counties. 2,189 transports (97%) were from lower-to higher-level facilities. Among these, 2,037 (93%) were to Level III tertiary care centers. 567 transports (25.2%) bypassed a closer facility offering an equivalent or higher level of care than its destination facility. Health system affiliation was associated with bypassing transport, indicating potential organizational rather than purely geographic drivers of transport decisions. Three "communities" of obstetric facilities largely shaped by geographic proximity were identified. ConclusionsAlthough Iowa does not have formal perinatal regions, patterns of maternal transport are mostly in line with three de facto regions. Some potential inefficiencies were identified, such as obstetric facilities transporting to a farther facility when a closer facility offered the same level of care or higher. These findings may help identify opportunities to enhance care coordination among obstetric facilities, optimize maternal transport networks, and improve regionalization of maternal care.

Background: Dengue remains a major public health challenge in Thailand despite decades of vector control implementation. While mathematical models have explored dengue transmission dynamics, systematic evaluation of current control strategies under realistic operational conditions remains limited. Methods: We developed a temperature-dependent, multi-serotype dengue transmission model that explicitly incorporates three primary vector control strategies: reduction in mosquito biting rates through personal protection measures, further reduction in mosquito birth rates beyond current larval control efforts, and further increase in adult mosquito mortality beyond current adulticide application levels. Using Approximate Bayesian Computation with Sequential Monte Carlo (ABC-SMC), we fitted the model to dengue hemorrhagic fever (DHF) surveillance data from nine province-year combinations representing high (Rayong), moderate (Ratchaburi), and low (Phrae) transmission settings across three years (2006, 2015, and 2017). The model accounts for four dengue serotypes, temperature-dependent mosquito dynamics, and temporary cross-protective immunity between serotypes. Results: The model closely reproduced observed monthly DHF case counts across all nine province-year combinations. Estimated reporting proportions ranged from 1.4% to 16.7%, with the highest values occurring in high-transmission provinces during the 2015 outbreak year. When each strategy was independently intensified by 50% relative to fitted baseline levels, reducing mosquito biting rates and increasing adult mosquito mortality consistently produced greater reductions in transmission than reducing mosquito birth rates. In the highest-transmission scenario (Rayong, 2015), a 50% reduction in biting rate from the baseline level yielded a 96.4% reduction in cumulative infections (95% CrI: 95.4-97.3%), compared with 94.3% (95% CrI: 91.8-95.6%) for a 50% increase in adult mosquito mortality and 77.0% (95% CrI: 58.6-84.6%) for a 50% reduction in mosquito birth rate. Analysis of the time-varying reproduction number (R_t) confirmed that interventions targeting adult mosquito-human contact achieved the greatest sustained epidemic suppression, although the relative ranking between bite prevention and adulticide application varied by epidemiological setting. Conclusions: Under the uniform 50% intensification scenario tested, interventions that directly disrupt adult mosquito-human contact, whether through personal protection or adulticide application, substantially outperformed larval control in reducing dengue transmission across diverse Thai settings. These findings support prioritizing personal protection and adulticide application, while the generalizability of this ranking to other intensification levels and settings warrants further investigation.

Species distribution models (SDMs) are widely used to support risk assessment for invasive non-native plant species (INNPS), but their performance is constrained by the coverage of occurrence data. Combining occurrences from citizen science (CS) platforms with data from structured state agency (StAg) monitoring provides unique advantages, yet they are rarely integrated. Here, we systematically compare how CS, StAg, and combined (COM) occurrence data influence the inferred environmental niches, predictive performance, and spatial applicability of SDMs for three widespread INNPS (A. altissima, H. mantegazzianum, I. glandulifera) in central Germany. We quantified niche overlap between datasets using PCA and Schoeners D and applied a hierarchical SDM utilizing boosted regression trees, while the Area of Applicability (AOA) was assessed to identify monitoring gaps. CS data were strongly biased toward lower-elevation, urbanized environments, whereas StAg data captured higher-elevation, remote habitats, particularly along watercourses. Niche overlap reflected both invasion stage and habitat preferences: A. altissima, a species that is spreading, showed the lowest overlap. H. mantegazzianum, associated with linear habitats like watercourses and infrastructure, exhibited intermediate overlap, while I. glandulifera, a widespread species, displayed the highest overlap. Overall, combined models achieved the highest predictive performance (AUC: 0.85, TSS: 0.58), reduced uncertainty along environmental gradients and produced more ecologically plausible suitability patterns. AOA analysis revealed high applicability ([&ge;]59%) across data sources and species, with COM models consistently reducing extrapolation uncertainty. Our findings highlight that integrating CS and StAg data reduces spatial biases and enhances SDM robustness, which is vital to improve INNPS risk assessments and management. HighlightsO_LICitizen science and state agency data capture distinct environmental spaces. C_LIO_LIOverlap between data sources is related to invasion stage and habitat preference. C_LIO_LICombined data improves invasive species niche representation and model accuracy. C_LIO_LIAOA analysis reveals monitoring gaps, especially in remote and high-elevation areas. C_LI

International guidelines for low-frequency electromagnetic field exposure (LF EMF) are primarily intended to prevent substantiated adverse effects. In the frameworks, limits on internal electric fields are linked to external exposure levels through computational dosimetry. However, the relationship between internal electric fields and these adverse effects remains incompletely understood. In particular, current approaches often overlook the morphological complexity and diversity of cortical neurons, which may limit the realism of neuronal activation estimates used to support these assessments. This study evaluates LF EMF-induced neural activation using 25 morphologically realistic neuron models spanning all cortical layers, embedded within 11 detailed human head models. The internal electric fields were simulated for uniform magnetic field exposures (100 Hz-100 kHz) along the three anatomical directions, and excitation thresholds were computed using a multi-scale framework combining voxel-based dosimetry with biophysical neuron simulations. A real-world exposure scenario involving a child near an acousto-magnetic article-surveillance deactivator was also analyzed. Thresholds varied across cell type, morphology, cortical location, subject anatomy, frequency, and exposure direction, with L2/3 pyramidal, L4 basket, and L5 thick-tufted pyramidal cells showing the lowest thresholds. Despite this variability, all simulated thresholds were conservative with respect to the basic restrictions and dosimetric reference limits set by IEEE ICES and ICNIRP. The smallest margin occurred at 100 kHz, where the threshold remained a factor of 2.8 above the corresponding limit. These findings indicate that current LF EMF exposure limits remain conservative when evaluated using highly detailed, morphology-based CNS activation models.

Canadas ambition for mineral security and its responsibilities to protect at-risk species and uphold Indigenous rights clash in the case of the Grays Bay Road and Port (GBRP) in Nunavut, an infrastructure project intended to unlock critical mineral deposits. We compiled Indigenous and Western science through a density analysis of caribou harvesting data near the proposed project site. We identified three consistently used harvesting hotspots, with the most significant hotspot lying directly in the path of the proposed GBRP project. These results indicate that the GBRP project will have significant and unmitigable negative effects on caribou conservation, food security, and Inuit harvesting rights. Prime Minister Carney claims that middle power countries must act consistently in this era of geopolitical rupture; this commitment must transfer to natural resource development reviews so that decision-making may be consistent and rooted in cross-legislation responsibilities and values, including the land claims agreements between Indigenous groups and the Government of Canada.

Background: Nigeria bears one of the highest maternal mortality burdens globally, with skilled birth attendance (SBA) remaining critically low in many regions. Understanding the independent determinants of SBA is essential for designing targeted interventions. Methods: This cross sectional study analyzed 21,465 births from the 2018 Nigeria Demographic and Health Survey (NDHS), a nationally representative household survey using stratified two stage cluster sampling. SBA was defined as delivery attended by a doctor, nurse, midwife, or auxiliary midwife. Multivariable logistic regression was used to estimate adjusted odds ratios (aOR) with 95% confidence intervals for the associations between SBA and maternal education, household wealth, place of residence, geopolitical region, maternal age, parity, and antenatal care (ANC) utilization, after accounting for confounding. Results: The overall prevalence of SBA was 44.9%. In the fully adjusted model, higher education (aOR = 7.01, 95% CI: 5.68-8.67), richest wealth quintile (aOR = 6.27, 95% CI: 5.27-7.46), and attending [&ge;]4 ANC visits (aOR = 3.80, 95% CI: 3.51-4.11) were the strongest independent predictors of SBA. Regional inequalities were pronounced, with SBA prevalence ranging from 17.7% in the North West to 85.6% in the South West. Crude effect estimates for education and wealth were substantially attenuated after adjustment, indicating large confounding by correlated socioeconomic factors. Conclusions: Maternal education, household wealth, ANC utilization, and geopolitical region are independent determinants of SBA in Nigeria. Scaling up ANC programs represents the most immediately actionable intervention, while long term gains require investment in girls' education and wealth equity. Targeted strategies for the northern regions are urgently needed. Keywords: skilled birth attendance, maternal mortality, Nigeria, DHS, antenatal care, logistic regression, health equity

Background: The increasing availability of routinely collected health data offers new opportunities for population-level research, yet access to comprehensive, linked, and standardised datasets remains limited. We describe EST-Health-30, a large-scale, population-representative health data resource from Estonia. Methods: EST-Health-30 comprises a random 30% sample of the Estonian population (~500,000 individuals), with longitudinal data from 2012 to 2024 and annual updates planned through 2026. Individual-level records are linked across five nationwide databases, including electronic health records, health insurance claims, prescription data, cancer registry, and cause of death records. A privacy-preserving hashing approach ensures consistent cohort inclusion over time while maintaining pseudonymisation. All data are harmonised to the Observational Medical Outcomes Partnership (OMOP) Common Data Model (version 5.4) using international standard vocabularies. Data quality was assessed using established OMOP-based validation frameworks. Results: The dataset contains rich multimodal information on diagnoses, procedures, laboratory measurements, prescriptions, free-text clinical notes, healthcare utilisation, and costs, with high population coverage and longitudinal depth. Data quality assessment showed high completeness and consistency, with 99.2% of applicable checks passing. The age-sex distribution closely reflects the national population, supporting representativeness, though coverage is marginally below the target 30% (29.2%), primarily attributable to recent immigrants without health system contact. The dataset enables construction of detailed clinical cohorts, analysis of disease trajectories, and evaluation of healthcare utilisation and outcomes across the life course. Conclusions: EST-Health-30 is a comprehensive, standardised, and population-representative real-world data resource that supports epidemiological, clinical, and methodological research. Its alignment with the OMOP CDM facilitates reproducible analytics and participation in international federated research networks, while secure access infrastructure ensures compliance with data protection regulations.

Background Polycyclic aromatic hydrocarbons (PAHs) and volatile organic compounds (VOCs) are combustion-derived pollutants linked to cardiovascular disease. Prior NHANES analyses have evaluated these chemicals individually, failing to capture the correlated co-exposure structures that characterize real-world environmental burden, thereby underscoring the need for application. In this study, we applied an unsupervised machine learning pipeline to urinary biomarker data to identify multi-chemical exposure clusters and quantify their differential cardiovascular risk profiles in a nationally representative US sample. Methods We analyzed 2,979 participants from NHANES between 2017-2018, representing an estimated 36.8 million US adults after complex survey weighting. Twenty-five urinary biomarkers (6 PAH, 19 VOC metabolites) were log-transformed, imputed using Multivariate Imputation by Chained Equations (MICE), and standardized. Uniform Manifold Approximation and Projection (UMAP) was used for dimensionality reduction, followed by Gaussian Mixture Model (GMM) clustering. Survey-weighted prevalence estimates with 95% confidence intervals (CIs) were calculated for hypertension and high total cholesterol within each cluster. Weighted multivariable logistic regression was used to estimate odds ratios (OR) for hypertension, adjusting for age, sex, race/ethnicity, and income. Results Four exposure clusters were identified with a mean assignment probability of 0.948. The High combustion cluster (n=370; estimated 5.1 million US adults) exhibited the highest multi-chemical burden and a weighted hypertension prevalence of 39.3% (95% CI 37.2-41.4%), compared to 28.7% (95% CI 21.9-35.5%) in the Low exposure reference group. After demographic adjustment, High combustion cluster membership was independently associated with 38.4% higher odds of prevalent hypertension (OR 1.38). The prediction model achieved a cross-validated area under the receiver operating characteristic curve (AUC) of 0.849 (SD 0.017). Non-Hispanic Black participants constituted approximately 40% of the High combustion cluster, exceeding their representation in lower-risk clusters. Conclusions Multi-chemical exposome profiling identifies four cardiovascularly distinct subpopulations in the US adult population. Membership in the High combustion exposure cluster was associated with higher odds of prevalent hypertension and disproportionately affected Non-Hispanic Black participants. These findings support the use of multichemical approaches over single-pollutant analyses and highlight the relevance of environmental exposure patterns for making policy and targeted cardiovascular risk stratification.

Ticks impose major health and economic losses on the livestock sector of Pakistan, yet uncertainty-aware maps of tick burden remain scarce. We focused on the two most common disease transmitting tick species, Rhipicephalus microplus and Hyalomma anatolicum, to produce exposure-adjusted district-level abundance estimates and predictions for unsampled areas in Punjab and Khyber Pakhtunkhwa (KPK). We compiled heterogeneous tick count records and standardized them per 100,000 animals. District-level climate and physiographic covariates were summarized via principal components analysis. Bayesian spatial models were fit in R-INLA using Gaussian likelihoods and BYM2 over a hybrid adjacency matrix. Competing non-spatial and spatial models were compared, and the best model was used to generate posterior predictions and 95% credible intervals for unsampled districts. Spatial models outperformed non-spatial alternatives and calibrated well. Model robustness was confirmed through eight independent 80/20 train-test splits, showing strong generalization with consistent predictions across seeds. For unsampled areas, R. microplus exhibited a pronounced north-south gradient with high predicted means but wide intervals in the northern highlands, indicating information gaps. H. anatolicum predictions were highest and most precise in southern Punjab. Sensitivity analysis highlighted a dominant spatial component, with modest contributions from PC1 and PC2. The Bayesian spatial models using the Besag-York-Mollie framework delivered comparable, exposure-adjusted tick abundance maps while quantifying uncertainty to guide surveillance. Results suggest a need for immediate control in confirmed hotspots and recommend targeted field sampling in high-uncertainty districts. The workflow generalizes to other vectors, pathogens, and regions for evidence-based livestock health planning.

Background Pseudomonas aeruginosa is a major cause of healthcare-associated infections in paediatric settings, where its persistence in moist environments such as hospital water and wastewater systems poses a particular risk to neonates and immunocompromised children. Aim The aim of this study was to showcase the long-term survival and transmission of P. aeruginosa in a large tertiary children's hospital in England which is crucial to develop strategies for water-safe care. Methods Environmental P. aeruginosa isolates were collected from taps, sinks, showers, and baths in augmented care areas of a 330-bed tertiary children's hospital built to NHS water-safety standards. Clinical isolates were classified as invasive (blood, cerebrospinal fluid, and bronchoalveolar lavage) or non-invasive (respiratory, urine, ear, abdominal, and rectal surveillance). Variable number tandem repeat (VNTR) profiles and metadata were extracted from PDF reports, de-identified, deduplicated, and curated using Python and R. Findings This retrospective study analysed nine-locus VNTR profiles of 457 P. aeruginosa isolates submitted to the UK Health Security Agency from a large tertiary children's hospital, identifying 56 isolate clusters (each with [&ge;]2 isolates), of which 19 (34%) contained at least one invasive isolate. The most persistent cluster (Cluster 1, n=20) spanned from July 2016 to September 2024, containing environmental and clinical (invasive and non-invasive) isolates. Conclusion These findings demonstrate long-term persistence of certain genotypes and temporal overlap between environmental and clinical isolates, highlighting the difficulty in detecting and eradicating P. aeruginosa in hospital water and wastewater systems and reinforcing the need for continuous rigorous water system controls.